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ABSTRACT 

The Seyfert 1 galaxy Zw 229-015 is among the brightest active galaxies being monitored by the Kepler 
mission. In order to determine the black hole mass in Zw 229-015 from H/3 reverberation mapping, we have 
carried out nightly observations with the Kast Spectrograph at the Lick 3 m telescope during the dark runs 
from June through December 2010, obtaining 54 spectroscopic observations in total. We have also obtained 
nightly V-band imaging with the Katzman Automatic Imaging Telescope at Lick Observatory and with the 
0.9 m telescope at the Brigham Young University West Mountain Observatory over the same period. We detect 
strong variability in the source, which exhibited more than a factor of 2 change in broad H/? flux. From cross- 
correlation measurements, we find that the V{j3 light curve has a rest-frame lag of 3 .86^0 days with respect to 
the V-band continuum variations. We also measure reverberation lags for Ha and H7 and find an upper limit to 
the H(5 lag. Combining the H/3 lag measurement with a broad H/3 width of cr\\ae = 1590 ±47 km s~' measured 
from the root-mean-square variability spectrum, we obtain a virial estimate of Mbh = 1-00;!;q 24 ^ 
the black hole in Zw 229-015. As a Kepler target, Zw 229-015 will eventually have one of the highest-quality 
optical light curves ever measured for any active galaxy, and the black hole mass determined from reverberation 
mapping will serve as a benchmark for testing relationships between black hole mass and continuum variability 
characteristics in active galactic nuclei. 

Subject headings: galaxies; active — galaxies: individual (Zw 229-015) — galaxies; nuclei 
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1. INTRODUCTION 

The NASA Kepler mission dBorucki et al.1 12010al) is de- 
signed to search for exoplanets by detecting transit events, and 
early results have demonstrated the spectacular quality of Ke- 
pler's ultra-precise photometrv (e.g. . iBorucki et al.l2010bb . In 
addition to its primary goal of obtaining light curves of Galac- 
tic stars, Kepler is monitoring a small number of active galac- 
tic nuclei (AGNs) of various types, including Seyfert galaxies, 
quasars, and blazars. Kepler will provide the best optical light 
curves ever obtained of AGNs, vastly exceeding the quality of 
any previous observations in terms of both sampling cadence 
and photometric accuracy. This will enable unprecedented 
measurements of the characteristic properties of AGN optical 
variability. 

One of the key motivations for measurement of high- 
precision AGN light curves is to determine the power spec- 
trum of continuum fluctuations and to search for any charac- 
teristic timescales of variability (in the form of breaks or slope 
changes in the power spectrum of variations), which might 
be related to light-crossing, dynamical, or thermal timescales 
in the accretion disk. At X-ray energies, long-duration light 
curves have been used to search for breaks in the power spec - 
tral density (e.g.. lUttlev et al.ll2002l: iMarkowitz et al.ll2()03h . 
and the break timescales are observed to be dependent on 
both black hole mass and luminosity across a range extend- 
ing from X-ray binaries to quasars (McHardy et al. 2006|. 
There is accumulating evidence that the optical fluctuations 
of AGNs also contain characteristic time scales that may de- 
pend on black hole mass and luminosity (jCoUier & Peterson! 
l200ltlKellv etani2009l: [MacLeod et alJl2oTol) ! 

Ground-based optical light curves of AGNs are limited 
in quality by nightly and seasonal sampling, weather, and 
seeing-dependent errors, but Kepler can provide exquisite 
data that will for the first time enable detailed measurements 
of the optical variability properties of individual objects over 
timescales ranging from hours up to (potentially) a few years. 
Unobscured (Type 1) AGNs monitored by Kepler can provide 
powerful tests of proposed relationships between variability 
properties, black hole mass, and luminosity that have been 
derived from larger samples of objects with lower-quality 
ground-based light curves. 

In order to fully exploit the results of these Kepler AGN ob- 
servations, independent measurements of black hole masses 
in Kepler-field AGNs are needed. Here, we present the results 
of a six-month, ground-based, optical reverberation-mapping 
program for the Kepler target Zw 229-015 (at z = 0.0275; 
iFalco et al.l 19991) . with the goal of determining the black hole 
mass in this galaxy from its H/3 reverberation lag. 

This galaxy (also known as CGCG 229-015 or 2MASX 
J19052592H-4227398) was first identified as a Type 1 Seyfert 
bv lProustl(ll990h . Since then, no follow-up observations of 
this AGN have been reported in the literature at all, aside 
from its inclusion i n ROSAT All-Sky Survey source catalogs 
dVoges et al.llI999l: Izimmermann et all 120011) . Nevertheless, 
Zw 229-015 is one of the nearest and brightest AGNs in the 
Kepler field, and it was selected as a monitoring target for the 
first two cycles of Kepler observations. To assess its suitabil- 
ity for ground-based monitoring, we obtained a test spectro- 
scopic exposure at the Lick Observatory 3 m Shane telescope 
in April 2010. This observation reconfirmed the Seyfert 1 
classification, and we then selected it for the reverberation 
mapping program reported here. In this paper, we describe 
our imaging and spectroscopic observations of Zw 229-015, 
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Fig. 1 . — V-band image centered on Zw 229-015, created by stacking the 18 
best-seeing frames (all with seeing better than 1 " 5) from the West Mountain 
Observatory. North is up and east to the left, and the size of the displayed 
region is 250" X 250" . Photometric comparison stars are circled and labeled, 
and the star that fell in the spectroscopic aperture is marked with a squai'e. 

the measurements of the continuum and emission-line light 
curves and the broad emission-line reverberation lag, and an 
estimate of the black hole mass based on the Hf3 lag and line 
width. 

2. IMAGING OBSERVATIONS 

Reverberation measurements require a light curve of the 
AGN continuum, and we obtained V-band photometry of Zw 
229-015 following t he same strategy used in recent rev erbera- 
tion programs (e.g., IWalsh et al.ll2009t iDenney et all 12010). 
Although ground-based measurements cannot approach the 
accuracy or sampling cadence of the Kepler data, they al- 
lowed us to monitor the target's variability continuously dur- 
ing the course of our observing program. Furthermore, the 
V-band photometry gives a nearly line-free measurement of 
the galaxy's continuum flux, while the broad Kepler spectral 
response includes the Ha and H/? lines from the AGN. 

We obtained regularly scheduled V-band imaging at the 
0.76 m Katzman Automatic Imaging Telescope (KAIT) at 
Lick Observatory, and at the 0.9 m telescope at the Brigham 
Young University West Mountain Observatory (WMO). At 
both telescopes, Zw 229-015 was observed each night dur- 
ing the monitoring period, except when precluded by weather 
or instrument problems. We also obtained a small number 
of images at other facilities, as described below. Figure [T] il- 
lustrates the field surrounding Zw 229-015. The galaxy has a 
barred spiral morphology, with a Hubble type of roughly SBa. 

2.1. KAIT 

The KAIT 0.76 m robotic telescope (Filipp enko et alj2001h 
uses a 512 X 512 pixel SITe CCD camera with a scale of 0."8 
pixel"'. A single 300 s exposure of Zw 229-015 was obtained 
on each night, weather permitting, from June 3 through De- 
cember 13 (UT dates are used throughout this paper). During 
the monitoring period, KAIT had two month-long shutdowns 
for repairs, first in August and again from mid-October to 
mid-November The median seeing for all KAIT exposures 
was 2." 15. The automatic KAIT data-processing pipeline in- 
cludes bias subtraction and flatfielding. 
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TABLE 1 
Photometric Comparison Stars 



Star 


a 


5 


V 




(h m s) 


(o / //) 


(mag) 


1 


19:05:19.49 


42:28:57.7 


14.983 ±0.006 


2 


19:05:35.95 


42:27:16.2 


16.575 ±0.014 


3 


19:05:33.71 


42:28:50.5 


15.954 ±0.011 


4 


19:05:19.98 


42:27:37.9 


16.090 ±0.016 


5 


19:05:27.28 


42:26:48.8 


16.000 ±0.015 


6 


19:05:30.05 


42:26:21.4 


16.198±0.016 



Note. — Coordinates are J2000. and are based 
on an astrometric solution obtained using the astrom- 
etiy.net package (Lang et al. 2010). The quoted mag- 
nitude uncertainties are calculated as the standard de- 
viation of all measurements from KAIT and WMO for 
the complete light curve of each comparison star 

2.2. West Mountain Observatory 

The 0.9 m telescope at the West Mountain Observatory is 
one of five small research telescopes operated by Brigham 
Young University. It was installed and tested in late 2009 and 
began working in a partially queue-scheduled mode during 
the 2010 observing season. The observations secured for this 
investigation were scheduled as part of an external guest in- 
vestigator program. 

From June 5 until September 1, observations were made us- 
ing a 3056 X 3056 pixel KAF 09000 CCD camera with a scale 
of 0."5 pixel"'. After September 2, observations were made 
using a 2048 x 2048 pixel Fairchild 3041-UV CCD camera 
with a scale of 0."6 pixel"'. The typical observing proce- 
dure was to secure two exposures of Zw 229-015 while it was 
close to meridian transit on each night that conditions permit- 
ted telescope operations. Exposures were 300 s from the start 
of monitoring until September 7, and were decreased to 200 s 
thereafter because of the higher sensitivity of the Fairchild 
CCD. Observations continued until December 13, but only a 
few exposures were obtained after November 5 due to poor 
weather Standard-star fields from Landolt (1992) were ob- 
served on a few photometric nights. The median seeing for the 
WMO exposures was 1 ."89. The standard processing pipeline 
for WMO data includes bias subtraction and flat-fielding. 

2.3. The Lick Nickel Telescope 

The Lick 1 m Nickel reflector was used to observe Zw 229- 
015 during a few nights while KAIT was shut down for re- 
pairs. The Nickel Direct Imaging Camera contains a thinned, 
Loral 2048 x 2048 pixel CCD with a scale of 0"184 pixel"', 
and the observations were obtained using 2 x 2 binn i ng for 
the CCD readout. On one photometric night, 'Landoltl (119921) 
standard-star fields were also observed. Overscan subtraction 
and flat-fielding using twilight sky flats were done following 
standard procedures in IRAFPI 

2.4. Maidanak Observatory 

Zw 229-015 was observed on five nights in August and 
September using the Seoul Nationa l University 4096 x 4096 
pixel camera (SNUCAM; llm et al.i,2010 ) and V-band filter 
at the Maidanak Observatory 1 .5 m telescope in Uzbekistan 
(Fairchild CCD; 0."266 pixel"'). The processing pipeline 

" IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 



for SNUCAM data included overscan subtraction and flat- 
fielding. 

3. PHOTOMETRIC REDUCTIONS AND MEASUREMENTS 

Our light-curve measurements followed the same proce- 
dures used by Walsh et al. (2009). First, each image was 
cleaned of cosmic-r ay hits using the LA-COSMIC procedure 
(I van Dokkuml200i1) implemented in IRAF. Aperture photom- 
etry on Zw 229-015 and several comparison stars was per- 
formed using the IRAF task phot. Ultimately, six compari- 
son stars were selected; these are illustrated in Figure [T] The 
photometric measurements were carried out with a range of 
aperture sizes. We found that the 4" aperture radius gave the 
cleanest light curves, and our final photometry is based on 
this aperture size. The background sky annulus had inner and 
outer radii of 8" and 14."4. When more than one exposure was 
taken per night at a given telescope, the measurements were 
averaged together to give a single data point for that night. 

Photometric calibration was accomphsh ed by carrying out 
aperture photometry on the iLandoltl (11992b standard stars ob - 
served at WMO and the Nickel telescope on nights that were 
judged to be photometric by the observers. We used only 
V-band exposures and did not attempt to derive any color- 
dependent terms for the transformation between instrumental 
and calibrated V magnitudes. We found excellent agreement 
between the photometric calibrations for the night of Septem- 
ber 26 at WMO and the night of August 12 at the Nickel tele- 
scope, and took the average of the two calibrations to obtain 
our final photometric scale. Table [T] lists the positions and 
magnitudes of the six comparison stars. 

The light curves from each telescope were separately nor- 
malized to V-band (Vega-system) magnitudes based on the 
comparison-star photometry. To check for small systematic 
differences between the photometric scales of the four tele- 
scopes, which could be present because color terms were 
not measured, we determined the average offset between the 
KAIT and WMO light curves for the 53 nights when the ob- 
ject was observed at both telescopes. The weighted average 
offset was 0.010 ±0.001 mag, and we applied this shift of 
0.01 mag to the KAIT light curve as an empirical correc- 
tion to bring it into agreement with the WMO photometric 
scale. Similar offsets of 0.010 and 0.020 mag were applied 
to the Nickel and Maidanak data points, respectively, in order 
to bring them into agreement with the combined KAIT and 
WMO light curve, based on the weighted average shift de- 
termined from nights when the galaxy was observed at more 
than one telescope. 

The final light-curve data for Zw 229-015 are listed in Ta- 
ble|2] Some observations taken in cloudy conditions had very 
low signal-to-noise ratio (S/N), and we excluded any mea- 
surement having a formal uncertainty greater than 0.03 mag. 
We also discarded measurements from a few images that suf- 
fered from very poor flat-fielding or other problems. The re- 
sulting light curve is illustrated in Figure |2] 

The photometric uncertainties listed in Table |2] are the val- 
ues returned by the photometry routine, which depend on 
photon-counting statistics for the object and background as 
well as on readout noise. These are likely to be underes- 
timates of the true uncertainties, since additional errors can 
result from flat-fielding irregularities, point-spread function 
variations, and other issues. One way to estimate the actual 
photometric uncertainties is to examine the magnitude fluc- 
tuations in the comparison-star light curves, which are pre- 
sumed to be intrinsically nonvariable. Comparison stars 3, 4, 
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TABLE 2 
Photometric Observations and 
Measurements 



HJD 


Telescope 


V 


-2450000 




(mag) 


5351.9785 


K 


15.989 ±0.008 


5352.8036 


W 


15.980±0.003 


5352.9688 


K 


15.976 ±0.007 


5353.9834 


K 


15.976 ±0.009 


5354.9619 


K 


15.957 ±0.008 



Note. — Telescopes are listed as follows: 
K = KAIT, W = WMO, N = Lick Nickel, M 
= Maidanak. This table is presented in its en- 
tirety in the electronic edition of the journal. 
A portion is shown here for guidance regard- 
ing its form and content. 

and 5 have V magnitudes very similar to that of Zw 229-015. 
For each of these stars we examined the scatter of measure- 
ments about the mean magnitude, and found standard devi- 
ations of 0.011-0.016 mag (see Table [T]!, while the median 
value of the formal photometric error on a single measurement 
for these stars is 0.006 mag. This discrepancy shows that the 
photon-counting statistics are usually not the dominant con- 
tributor to the error budget, except when the images have very 
low S/N. We address this issue by adding an additional error 
term in quadrature to the photometric uncertainties prior to 
carrying out the cross-correlation measurements, as described 
ing6] 

Our photometric aperture includes a substantial portion of 
the host galaxy, including essentially the entire bulge and 
much of the bar and disk. We do not attempt to subtract 
off this host-galaxy starlight contribution to the light-curve 
points, because our ground-based imaging does not have suf- 
ficiently high resolution to decompose the AGN point source 
and bulge into separate components. Removing this constant 
contribution of host-galaxy starlight would not affect the mea- 
surement of the reverberation lag. 

A condensed version of the V-band light curve was pre- 
pared for the cross-correlation measurements by taking a 
weighted average of any photometric data taken within 5 
hours of each other This produces a light curve with at most 
one data point per night. The condensed light curve was then 
converted from magnitudes to linear flux units. 

4. SPECTROSCOPIC OBSERVATIONS 

The spectroscopic monitoring program at the Lick 3 m 
Shane telescope was carried out by several independent ob- 
serving teams. Unlike other recent reverberation datasets that 
were obtained via long-duration, dedicated campaigns of clas- 
sically scheduled observing nights, we chose to obtain synop- 
tic monitoring of this single target by enlisting the entire com- 
munity of observers using the Kast spectrograph at Lick dur- 
ing Summer and Fall 2010. The Kast spectrograph is mounted 
at the Cassegrain focus of the 3 m telescope for about half 
of each lunation, and the spectrograph is used for a variety 
of programs including supernova and transient follow-up ob- 
servations, stellar spectroscopy, and AGN surveys. During 
each Kast observing night, the regularly scheduled observers 
took a 1200 s exposure of Zw 229-015 to contribute to this 
monitoring program. This strategy enabled us to collect a 
spectroscopic dataset with almost nightly sampling during 
each dark run, except for weather-related losses. The spec- 
troscopic monitoring began with one observation in May, and 
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Fig. 2.— y-band light curve of Zw 229-015. 

then continued with nightly observations during each dark run 
(weather permitting) from June through mid-December From 
June through September, nearly every night was clear enough 
to observe, while October through December had more sub- 
stantial weather losses, as is typical for Mt. Hamilton (the site 
of Lick Observatory). Table |3]gives a log of the Kast obser- 
vations for the 54 nights when useful data were obtained. We 
have omitted data from two nights when Zw 229-015 was ob- 
served through very thick clouds, resulting in a S/N that was 
too low to be useful for any measurements. 

Each observing team used a different instrumental setup 
for the spectrograph, designed for their primary science pro- 
grams. On the blue side of the spectrograph, all observers 
used the 600-line grism, which gives a dispersion close to 
1.0 A pixel"'. The spectral range falling on the blue-side de- 
tector is controlled by adjusting the CCD position via a man- 
ually movable stage, resulting in a fixed wavelength coverage 
for the entire night's observations. Most observing teams set 
the CCD to a position that covered 3440-5520 A, placing 
the H/3 and [O III] AA4959, 5007 lines on the blue CCD. How- 
ever, on several nights a bluer setup covering ^ 3 140-5 150 A 
or ^ 3125-5135 A was used, and with this setup the [O III] 
A5007 line fell off the edge of the blue CCD, while the A4959 
line was still present. 

For the red side of the spectrograph, several gratings are 
available and different gratings were installed on different 
nights. Since we could not ensure the same degree of con- 
sistency in the red-side observations as was obtained on the 
blue side, the red-side observations were a lower priority for 
this project. On some nights Zw 299-015 was not observed 
at all on the red side. When it was observed, gratings of 300, 
600, and 1200 lines mm"' were used for the red-side observa- 
tions, with dispersions ranging from 1.2 to 4.6 A pixel"'. In 
most cases, the observers used a d55 dichroic to separate the 
blue and red beams of the spectrograph, while on a few nights 
a mirror was used that directed all of the light to the blue side. 
When the 300 or 600 lines mm"' gratings were used, there 
was usually a small region of wavelength overlap between the 
blue and red spectra. As a result of the heterogeneous nature 
of the red-side data, this paper focuses primarily on the anal- 
ysis of the blue-side observations, which were all taken with 
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the same grism and homogeneously observed except for the 
variations in spectral range described above. 

With a single exception (the night of 2010 August 16), all 
observations of Zw 229-015 were obtained with the spectro- 
graph slit at a position angle of PA = 56. °5. This angle was 
chosen so as to place a bright foreground star (marked with a 
square in Figure [T) in the slit, in order to aid in centering the 
primary target. A 4" slit width was used for all exposures of 
Zw 229-015 and for the standard stai-BDH-28°4211. A 1200 s 
exposure of Zw 229-015 and a 60 s exposure of BDh-28°421 1 
were taken on most nights. On a few nights, the observers 
adjusted the exposure time due to weather constraints or took 
more than one exposure of the AGN. The airmass was smaller 
than 1 .2 for 72% of the spectroscopic observations of Zw 229- 
015, and smaller than 1.5 for all but the final three exposures, 
when the object was observed at airmasses between 1.5 and 
2.0. At the end of the observing season in December, when 
the target was observed at relatively high airmass, the fixed 
slit PA of 56.°5 was wit hin 10-20° of the optimal parallactic 
angle dFilippenkoll 1 9821) for the exposures, resulting in mini- 
mal differential light losses. 

Table [3] also lists the "image quality," measured as the full 
width at half-maximum intensity (FWHM) of the spatial pro- 
file of the Galactic foreground star in the spectroscopic ex- 
posure. This image quality is primarily a measure of atmo- 
spheric seeing, but on some nights poor focus or occasional 
guiding errors added significantly to the spatial profile width. 
On clear nights, the blue-side observations typically resulted 
in a S/N per pixel of ~ 40^5 in the extracted spectra at 
4700 A. 

Dome-flat exposures and observations of arc-line lamps 
were taken during each afternoon. Dome flats were usually 
observed through a 2" slit, since that slit width was used for 
the observers' primary programs. The difference in slit width 
between the AGN exposures and dome-flat exposures did not 
have any significant impact on the quality of the blue-side flat- 
tening, but on the red side it left strong fringing residuals at 
wavelengths longer than ^ 7000 A. 

5. SPECTROSCOPIC REDUCTIONS AND MEASUREMENTS 

Reduction of the Kast data followed standard procedures, 
including overscan subtraction, flat-fielding, unweighted 
(nonoptimal) extractions, and wavelength calibration using 
the line-lamp exposures. An extraction width of 6."5 was 
used. Flux calibration and a linear shift to the wavelength 
scale based on the wavelengths of night-sky emission lines 
were then applied to the spectra. The final blue-side spectra 
were binned to a uniform linear scale of 1 A pixel"'. Cosmic- 
ray residuals were removed from the reduced spectra by inter- 
polation over the affected pixels. For nights when more than 
one exposure was taken, the multiple exposures were com- 
bined into a single reduced spectrum. Error spectra were also 
extracted and propagated through the full sequence of reduc- 
tions. Figure |3] illustrates the reduced spectrum from a single 
observation. 

We measured the centroid wavelength of the [O III] A5007 
line in each spectrum where the line was available. The me- 
dian result is Aobs = 5143. 37± 1.17 A, corresponding to an 
emission-line re dshift of z = 0.02 73 ± 0.0002, consistent with 
the value from .Falco et al.l (119991) . The fairly large dispersion 
in the [O III] wavelength measurements is most likely a result 
of variations in the centering of the AGN in the wide 4" slit 
in each observation. 
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Fig. 3. — Combined Kast blue-side and red-side spectram of Zw 229-015, 
from 2010-06-13. On the red side, a 1200-lines mm"' grating was used, 
giving a dispersion of 1.2 A pixel"' . 

In order to intercalibrate the flux scales of the spectra and 
align them to a consistent waveleng th scale, we employed the 
scaling procedure described by v an Groningen & WandersI 
(1992), which is based on the assumption that the narrow 
[O III] emission-line flux is intrinsically constant over the 
timescale of the monitoring program. This method applies 
small wavelength shifts and a multiplicative flux scaling to 
each spectrum, as well as a convolution with a Gaussian 
broadening kernel, and finds the values of the shift, kernel 
width, and flux scaling factor that minimize the residual dif- 
ferences between the individual spectrum and a reference 
spectrum over a small wavelength region containing a nar- 
row emission line and some surrounding continuum. Accord- 
ing to lvan Groningen & WandersI (119921) . their scaling method 
can yield accuracies of 5% or better, but higher accuracy 
(r^ 1-2%) can be obtained with the be st-quality datasets (e.g., 
iBentz et alj|2009l: iDenney et alJl201 0'). The flux scale of the 
reference spectrum was normalized by scaling it to match the 
mean [O III] A5007 flux measured from eight nights that were 
considered to be photometric by the observers. For those 
eight nights, we find a mean observed [O III] A5007 flux of 
3.96 X lO"''* erg cm"^ s"\ with a standard deviation of 4%, 
as measured from the flux-calibrated spectra prior to applying 
the scaling procedure.^ 

When applying the Ivan Groningen & WandersI (119921) pro- 
cedure to the individual spectra, we used the [O III] A4959 
emission line as the flux scaling reference line, instead of the 
A5007 line. The reason for this modification was that 10 of 
the spectroscopic observations used the blue-side CCD po- 
sition that cut off the A5007 emission line. The wavelength 
range in common to all of the blue-side exposures extended 
to 5130 A, including the A4959 emission line and about 20 A 
redward of the line. By using the A4959 line as the flux refer- 
ence, all of the spectra could be treated identically in the flux 
scaling procedure. 

To examine the accuracy of the spectral scaling, we mea- 
sured the flux in the A4959 emission line from each spectrum 
after the scaling was applied, along with its uncertainty from 
the propagated error spectrum. We then determined the nor- 
malized excess variance al of the A4959 light curve, follow- 
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Fig. 4. — Mean and RMS spectra. The plot illustrates the wavelength range 
in common to all of the blue-side spectra over the course of the campaign. 
Narrow emission-line residuals in the RMS spectrum become stronger toward 
the blue end due to errors in wavelength calibration and differences in spectral 
focus between nights, since the spectral scaHng optimizes the alignment and 
resolution of the spectra based on the A4959 emission line at the red end of 
the blue- side spectra. 

ing the definition used by lNandra et alJ (Il997h : 



(1) 



where is the total number of observations, /i is the mean 
flux, and X, and cr, are the individual flux values and their un- 
certainties. We find CTx = 0.024, indicating that the residual 
errors in flux scaling contribute a spurious scatter to the light 
curves at the ^2% level. These residuals are primarily due to 
seeing variations (particularly from the worst-seeing nights), 
miscentering of the target in the slit, and night-to-night varia- 
tions in the telescope and spectrograph focus. Since Zw 229- 
015 exhibited dramatic, high-amplitude variability during the 
monitoring period, this residual scatter only has a small effect 
on the H/3 light curve. 

Figure |4] illustrates the mean and root-mean-square (RMS) 
spectra based on the full dataset, constructed after applying 
the spectral scaling procedure to the data. The RMS spec- 
trum essentially gives the standard deviation of flux values at 
a given wavelength, and illustrates the relative degree of vari- 
ability of different wavelength regions of the spectrum. The 
broad H/3, H7, and lines clearly stand out in the RMS 
spectrum, and the variable continuum rises strongly toward 
the blue end of the observed wavelength range. In addition, 
a broad He II A4686 line is evident in the RMS spectrum, 
although it is too weak to be visible in the mean spectrum, 
where only the narrow component is clearly detected. 

The flux scaling of the red-side data is more difficult to es- 
tablish accurately. On the red side, t here are no strong and iso- 
lated n arrow emission lines, so the van Groningen & WandersI 
(11992b scaling routine cannot be applied directly. Addition- 
ally, on some nights the standard star was not observed si- 
multaneously on the blue and red sides, resulting in flux scal- 
ing offsets when conditions were not photometric or the see- 
ing was poor In order to place the red-side data on a flux 
scale that was approximately consistent with the blue-side 



data, we first scaled each red spectrum by the scale factor 
used for the corresponding blue spectrum. This worked well 
for most nights, but left obvious offsets between the blue-side 
and red-side flux scales for a few nights. For those nights, 
an additional scaling factor was applied to the red-side data 
in order to bring the spectrum into agreement with the cor- 
responding blue-side spectrum over the region of wavelength 
overlap or near-overlap. (We note that in a dedicated moni- 
toring campaign with homogeneous observations, these flux 
scaling problems can be mitigated by using a spectrograph 
setup with substantial wavelength overlap between the red and 
blue sides, and by observing standard stars simultaneously on 
both arms of the spectrograph.) 

The H/3 light curve was measured by first subtracting a lo- 
cal, linear continuum underlying the line. The continuum fit- 
ting windows were set to 4860^910 A and 5110-5130 A. 
The H/3 flux was then measured by direct integration of the 
spectrum over 4930-5050 A. We similarly measured the light 
curves of H7 and H(5. The emission-line summation windows 
for H7 and H5 were 4420-4500 A and 4175-4245 A, respec- 
tively, while the continuum windows were 4300^400 A and 
4550^600 A for H7, and 4100-^150 A and 4275^325 A 
for H(5. On the red side, the emission-line window for Ha 
was 6600-6825 A, and the continuum windows were 6500- 
6560 A and 6940-7020 A. AH of these emission-fine mea- 
surements include contributions from the narrow components 
of the Balmer lines, which are presumed to be nonvariable 
over the timescale of the monitoring project. In the mean 
spectrum, the narrow component of H/3 contributes ^ 10% 
of the total H/3 flux. In addition, the H7 integration window 
includes the [O III] A4363 line, and the Ha window includes 
[N II] AA6548, 6583. We do not attempt to decompose these 
emission blends into separate components. 

Figure|5]displays the Balmer-line light curves. The H/3 light 
curve is the cleanest of the four, because of the proximity of 
the H/3 line to the [O III] line used as the flux scaling refer- 
ence. The Ha light curve is substantially noisier than that of 
H/3, despite the higher flux of Ha, due to the issues in match- 
ing the flux scales of the red-side and blue-side spectra. Nev- 
ertheless, all of the lines clearly show the same pattern of flux 
variations throughout the monitoring period. 

The variability amplitude of the light cur ves can be quanti- 
fied using the statistics j^max and /y^r ( e-g-, iKaspi et alJl2000l : 
iPeterson et aLll2004t iBentz et al.ll2009l) . where iJ^ax is the ra- 
tio of maximum to minimum flux observed, and ivai is calcu- 
lated according to 



(2) 



where and (5^) are the variance and the mean-square mea- 
surement uncertainty in the fluxes, and (/) is the average flux 
in the light curve. For the H/3 light curve of Zw 229-015, 
we find /?max = 2.68 ± 0.04 and /Var = 0.23. By both of these 
measures, Zw 229-015 displayed strong variability, exceeding 
the values of Rmax ™d /\,ar observed for the entire sample of 
13 AGNs observed in the 2008 Lick AGN Monitoring Project 
during a shorter 2-month campaign (Bentz et al. 2009). 

Figure |5] also displays the light curve of the continuum re- 
gion 3550-3800 A. (We refer to this spectroscopic f/-band 
continuum measurement as the Us light curve.) This is the re- 
gion most dominated by the AGN continuum in the blue-side 
spectra. The spectroscopic continuum fight curve is consid- 



REVERBERATION IN ZW 229-015 



7 



15.9 
(0 16 

B 

16.1 



16.2 
1.2 
1 

^ 0.8 

0.6 ^ 3550-3800 A 



200 
180 
160 
7 140 
^ 120 



o 40 
30 



20 

20 
15 
10 

10 
8 
6 
4 



— 1 1 1 1 1 r- 



,1 . 



I'm 



' 



I,, II. , II I' 

V .,' ' I'V" ' ' I 

V I ,1 I I i| ' 
\ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 H 



I','' 



I 



I 'I 



I |l| ' 



11, I 



H 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 H 



Ha 



H h 



H 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1- 



H h 



H 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1- 



H 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 h 



" 



; H6 

_l I I I I I I I I L_ 



5350 5400 5450 5500 

HJD - 2450000 



5550 



Fig. 5. — Light curves of the V-band continuum, the spectroscopic continuum over 3550-3800 A (in units of 10 erg cm ^ s ' A '), and the Ha, H/3, H7, 
and H<5 emission lines. The displayed V-band light curve is the "condensed" version in which measurements separated by less than 5 hr have been averaged 
together giving a single data point per night. The 3550-3800 A continuum light-curve points are plotted with error bars of ±3% as a rough approximation of the 
unceitainties from the spectral scaling procedure and slit losses, since these en'ors dominate over photon-counting uncertainties. 



erably noisier than the V-band light curve, as expected due 
to the combination of seeing-dependent slit losses, miscenter- 
ing of the AGN in the slit, and the nonparallactic slit orien- 
tation. Nevertheless, the measurements illustrate that the U - 
band continuum generally follows the same pattern of varia- 
tions as the V-band light curve, with a factor of ~ 2 variability 
amplitude during the monitoring period. 

6. REVERBERATION MEASUREMENTS 

6. 1. The Hj3 Emission Line 

Although it is difficult to see in Figure|5] the continuum and 
emission-line light curves do exhibit a small temporal offset. 
To demonstrate this more clearly. Figure |6]illustrates the cen- 
tral portion of the V-band continuum light curve (in flux units) 
and the H/3 light curve, after rescaling to approximately the 



same flux level. Aligning the two light curves by eye, it is 
evident that the H/3 variations lag the continuum by ^ 4 days. 

To determine the reverberation lag quantitatively, we em- 
ployed the interpolation cross-correlation function (ICF) 
method (Gaskell & Peterson 1987), f ollowing the prescrip- 
tions described by White & P eterson! (Il994l) . iPeterson et al.l 
(.2004.) . and .Bentz et al. (.2009). The ICF method deals with 
unevenly sampled time series by interpolating linearly be- 
tween data points to obtain temporally matched time series for 
the continuum and emission-line light curves. For a given trial 
value of the lag time r, the continuum light curve is shifted 
forward in time by r days, and the cross-correlation is com- 
puted for the region of temporal overlap between the two light 
curves. The value of the cross-correlation function (CCF) for 
a given lag is computed twice: first by interpolating the shifted 
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TABLE 3 

Spectroscopic Observations and Measurements 



UT Date 


HJD 


FWHM 


S/N 


ma) 


/(H/3) 


/(H7) 


/(H<5) 




-2450000 


(arcsec) 






(10-'5 erg 


cm ^ s ' ) 




2010-05-19 


5335.995 


2.2 


43 




42.39 ±0.24 


18.27 ±0.20 


8.54±0.18 


2010-06-08 


5355.978 


2.1 


45 


184.52 ±0.62 


44.69 ±0.25 


20.44 ±0.20 


10.69 ±0.19 


2010-06-11 


5358.954 


2.9 


37 


190.32 ±0.30 


43.19±0.32 


19.56 ±0.26 


10.56 ±0.25 


2010-06-13 


5360.985 


2.5 


40 


190.51 ±0.23 


40.25 ±0.27 


17.76 ±0.22 


8.12±0.21 


2010-06-15 


5362.928 


2.7 


40 


181.04±0.33 


38.19±0.29 


17.04 ±0.24 


9.11 ±0.22 



Note. — The listed FWHM image quality is measured from the spatial profile of the bright star that 
fell on the spectrograph slit in each observation. The listed S/N is the signal-to-noise ratio per pixel in the 
reduced blue-side spectrum at A = 4700 A. Measured fluxes include the contribution of blended narrow-line 
components. This table is presented in its entirety in the electronic edition of the journal. A portion is shown 
here for guidance regarding its form and content. 



continuum light curve to the time steps of the emission-line 
light curve, and then by interpolating the emission-line light 
curve to the time steps of the shifted continuum light curve. 
The two interpolated versions of the CCF are then averaged 
together to give the final CCF. We computed the CCF for lag 
values ranging from -20 to 30 days, in increments of 0.25 
days. 

Once the CCF is computed, the emission-line lag is deter- 
mined in two ways: the peak of the CCF defines Tpeak, while 
the centroid of the CCF for points above 80% of the peak 
value of the CCF gives Teen- The value of Tpeak and its un- 
certainty are quantized in units of the sampling interval used 
for the cross-correlation, in this case 0.25 days. Black hole 
masses from reverberation data are most often determined us- 
ing Teen (Peters on et al .1120041) . 

The final cross-correlation results and their uncertainties 
were determined by follo wing the Monte Carlo bootst r appin g 
pro cedure descri bed by P eterson et all (11998b . IWelshI (119991) . 
and lPeterson et al. (2004). In this method, a large number of 
modified realizations of the light curves are created and the 
cross-correlation function is computed for each realization of 
a pair of continuum and emission-line light curves. Each sim- 
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Fig. 6. — Central portion of the V-band continuum and H/3 light curves, 
illustrating the reverberation lag. The light curves have been scaled to the 
same overall flux level, and a constant offset was first subtracted from the 
continuum so that its variability amplitude would match that of the H/3 light 
curve. In the lower panel, the H/3 light curve is shifted by an offset of -4 
days. 



ulated light curve is generated by selecting n points randomly 
from the original light curve, where n is the number of points 
in the original light curve and the random selection allows a 
particular point to be chosen more than once. If a particular 
point is selected m times, then its uncertainty is reduced by 
a factor of m'/^. Then, the simulated light curve is deviated 
by adding random Gaussian noise based on the uncertainty 
in each data point. The emission-line and continuum light 
curves are each resampled 10"* times by this procedure, and 
for each pair of simulated light curves the cross-correlation 
function is determined. From these simulations, we build up 
distributions of values of Teen and Tpeak- The final adopted val- 
ues of Teen and Tpeak are defined to be the median values from 
each distribution. The Icr uncertainty ranges on Teen and Tpeak 
are given by the range of lag values around the median that 
exclude the highest 15.87% and lowest 15.87% of lag values 
in the distribution. 

When the continuum and emission-line light curves have 
very similar observational cadences, the two interpolated ver- 
sions of the CCF are typically very similar in shape. However, 
if the two light curves have very different temporal sampling, 
then the two CCFs can d iffer substantially. As described by 
IWTiite & Petersonl(ll994l) . the two interpolated versions of the 
CCF are usually averaged together to obtain the final result, 
unless there is a strong reason to prefer one interpolation over 
the other We tested whether this might affect the lag mea- 
surements by measuring Tpeak and Teen separately for the two 
interpolated versions of the CCF, and comparing them with 
the results obtained for the averaged CCF. Although the two 
interpolated CCFs differed somewhat in shape, the measured 
lags were very similar in both cases, and the derived values 
of Teen differed by only 0.18 days for the two separate CCFs. 
Given this close agreement, we chose to use the averaged CCF 
for our final results. 

As discussed by Welshl (119991) . the reliability of cross- 
correlation lag measurements can be improved by subtract- 
ing off any long-term secular variations from the light curves 
that occur on timescales much longer than the reverberation 
timescale. In practice, this is often accomplished by subtract- 
ing a simple linear fit from the light curves. We found that the 
measurements of Tpeak and Teen were virtually unchanged re- 
gardless of whether the light curves were linearly detrended, 
and detrending with a quadratic fit (in an attempt to remove 
the long-timescale curvature of the light curves) only mod- 
ified the measured lags by less than their Icr uncertainties. 
For our final measurements, we chose to use the linearly de- 
trended light curves. We also tested to see whether the cross- 
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Fig. 7. — Top: Autocorrelation function of the V-band light curve. Bottom: 
Cross-correlation function of the H/3 and V-band light curves. The dotted 
and dashed lines mark the values of t^-^ and Tccn, respectively. 

correlation measurements would be affected by removal of the 
last few spectroscopic data points, during the period of poor 
weather when the sampling was very sporadic, and found that 
the results for Tpeak and Teen again were consistent with our 
final adopted values within the Icr uncertainties. 

As described above in ^ the photometric errors listed 
in Table |2] are likely to underestimate the true uncertain- 
ties. To account approximately for additional sources of er- 
ror, we added 0.01 mag in quadrature to all of the V-band 
photometric errors prior to computing the cross-correlations. 
We found that this change had only a small effect on the 
cross-correlation results: the derived value of Teen was either 
3-85;!;ogQ or 3.97;!;Qg3 days using the original photometric er- 
rors or the expanded error bars, respectively. 

Our final measurement of the CCF gives observed-frame 
H/3 lags of Teen = ^ -91^11 days and Tpeak = 3.50;!];50 days, 
consistent with the general expectation of a ^ 4 day lag from 
inspection of the light curves (as in Figure |6]l. Converting to 
the rest frame of Zw 229-015, the lags are Teen = 3.86;!;q and 
T'peak = 3.41^0 71 days. All of the lag measurements are listed 
in Table |4] Figure |7] displays the autocorrelation function of 
the V-band light curve and the H/3 vs. V-band CCF. The peak 
value of the CCF is 0.86, indicating a fairly strong correlation 
between the continuum and emission-line variations. 

We also measured the H/3 lag relative to the f/j-band light 
curve and found observed-frame lags of Teen = 5.16^[ and 

^peak 

= 3.25^^0 75 days. Within the substantial uncertainties, 
these results are consistent with the lag measurements carried 
out against the V band. Cross-correlating the V-band light 
curve against the continuum, we find values of Tpeak and 
Teen that are consistent with zero within their 1 a uncertainties, 
so there is no significant detection of a lag between contin- 
uum bands. Higher-cadence observations would be needed 
in order to search for any cross-band continuum lags (e.g., 
ISergeev et al.ll2005l) . In addition to the featureless continuum 
of the AGN, the spectral region also includes emission 
from the "small blue bump," which contain s both B aimer con- 
tinuum and Fe II emission. For NGC 5548. lMaoz et alj (119931) 
showed that the small blue bump had a lag relative to the UV 
continuum that was similar to the lag of Ly a, consistent with 



the expectation that the small blue bump emission originates 
from the broad-line region (BLR). The nondetection of a lag 
between the Us and V bands in Zw 229-015 suggests that ei- 
ther the small blue bump makes only a small contribution to 
the Us band, or that the Balmer continuum has an extremely 
short lag time relative to the primary continuum. The latter 
possibility would be consistent with the decreasing lags mea- 
sured for higher-order Balmer lines (see § l6.3l below). 

There is a small contribution to the V-band flux from the 
broad H/3 line, and also an extremely small contribution from 
Ha at the red end of the Johnson V filter passband. This can 
potentially bias the lag measurement, by adding a component 
to the continuum light curve that has zero temporal shift with 
respect to the emission-line light curve. Consequently, the 
lag measured from the cross-correlation peak will tend to be 
lower than the true lag. In order to assess the magnitude of 
this effect, we carried out synthetic photometry on Zw 229- 
015 using a combined blue-side and red-side spectrum con- 
structed from a night with good overlap between the blue and 
red sides. Using the IRAF S YNPHOT package, we calculated 
the V-band magnitude of the spectrum twice: first on the orig- 
inal spectrum, and then on a modified version of the spectrum 
in which the H/3 and Ha emission lines were completely re- 
moved by manually interpolating over them. The two syn- 
thetic V-band measurements differed by only 0.02 mag, in- 
dicating that the broad Balmer lines make only a very small 
contribution (s^ 2%) to the broad-band V magnitude. 

To obtain a quantitative estimate of the magnitude of this 
bias on the lag measurement, we carried out a set of Monte 
Carlo simulations. We simulated a larg e number (10"*) of 
AGN light curves using the method of iTimmer & Koenigl 
dun, for a power-density spectrum of the form P(f) oc / 
normalized to a continuum RMS variability of 20% over a 
180-day period. A simple delta-function transfer function 
with a lag of 4 days was used to model the emission-line re- 
sponse to the continuum fluctuations. The "contaminated" 
V-band light curve was simulated by adding the emission- 
line light curve to the continuum light curve, weighted so 
that 2% of the total flux was from the emission-line contribu- 
tion. With each set of simulated light curves we measured the 
emission-line lag relative to the uncontaminated continuum 
light curve, and relative to the continuum light curve with the 
added emission-line contribution. 

The simulations confirm that the bias in the lag measure- 
ments is very small forZw 229-015. When the emission-line 
contribution to the continuum light curve is at the observed 
2% level as in Zw 229-015, the simulations show that the 
measurement of Teen is biased on average by only -0.07 days 
(with an RMS scatter of 0.05 days in the bias for the full en- 
semble of simulations), relative to the measurement obtained 
by cross-correlating against the pure continuum light curve. 
For Zw 229-015, the bias is much smaller than the overall un- 
certainty in the lag measurement, but it can become significant 
for higher levels of emission-line contamination. We carried 
out additional simulations, and found that if the emission-line 
contribution to the broad-band light curve increases to 10%, 
for example, the average bias in the Teen measurement is -0.3 
days (for the same model parameters described above). For 
future reverberation campaigns that use broad-band photome- 
try to measure the continuum light curves, it will be important 
to carry out simulations like these to assess the likely magni- 
tude of this effect on each individual object. The size of the 
bias would depend on both the emission-line equivalent width 
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TABLE 4 
Lag Measurements 



Measurement Tccn (days) Tp^ak (days) 



Ha vs. V 
vs. V 
H7 vs. V 
ns vs. V 
H/3 vs. Us 
V vs. Us 



5 22+0-'*3 

•^•^ -0.93 
5 16+1 -2' 

-1.76 



0, 



4 75+1.00 

3 5()+''-5<' 

2 SO+O '^ 
-0.25 

^- '-^-0.50 

-0 25+^-^5 



Note. — All lags are in the observed 
frame. Rest-frame lags can be obtained by di- 
viding by l+z, or 1.0275 in this case. Us de- 
notes the C/-band continuum measured from 
the blue-side spectra over 3550-3800 A. 

and on the target's redshift, since higher redshifts would move 
the broad H/3 line closer to the peak transmission of the V fil- 
ter Alternatively, another approach would be to subtract off 
the emission-line flux directly from the broad-band light curve 
before measuring the cross-correlation lag. 

6.2. Velocity-Resolved Measurements 

The behavior of the reverberation lag as a function of veloc- 
ity across broad emission lines is dependent on the kinematics 
of the BLR, and recent ground-based programs have success- 
fully obtained velocity-resolved measurements of reverbera- 
tion lags for H/3 and other lines in several AGNs (Kollatschnv 
[2003; Bentz eiaL 2009, 2 010a; Dennev et al. 2010). Our Zw 
229-015 dataset is not ideally suited for velocity-resolved re- 
verberation measurements, because the H/3 lag is fairly short 
and the monthly gaps in the spectroscopic light curve limit 
the measurement accuracy for lags of individual velocity seg- 
ments of the H/3 line. 

Nevertheless, we are still able to obtain some rudimentary 
information on the velocity dependence of the reverberation 
lag across the H/3 line profile. Light curves were measured 
for several independent velocity segments across the H/3 line, 
and the individual segment light curves were cross-correlated 
against the V-band continuum. The results are illustrated in 
Figure [8] Six velocity segments were used for the final mea- 
surements. In the line wings, where the variability amplitude 
is relatively low, the measurements have large uncertainties; 
the lags measured for the blue-wing and red-wing segments 
are consistent with zero within their la uncertainties. (We 
used a broader velocity segment for the blue wing because 
subdividing the blue wing into even smaller segments led to 
very poor cross-correlation results, with some subsegments 
having negative lags.) 

With these caveats in mind, the velocity-resolved behav- 
ior seen in Figure |8] shows marginal evidence for a slightly 
longer lag in the line core relative to the mean lag, while 
the lag in the line wings is shorter A symmetric blue-red 
response in the velocity-delay map, with shorter lags in the 
line wings, is qualitativel y consistent with Kepl erian motion 
of the BLR clouds (e.g.. IWelsh & HomelfT99lb . but due to 
the large uncertainties in the velocity-resolved measurements 
we refrain from drawing any specific conclusions about the 
BLR kinematics. The results do illustrate that Zw 229-015 is 
a promising candidate for further velocity-resolved reverbera- 
tion work. Given the strong and rapid variability in this object, 
a high-cadence monitoring campaign could obtain much bet- 
ter results, and a reconstruction of the two-dimensional trans- 
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Fig. 8. — Velocity-resolved reverberation results for H/9. The lower panel 
shows the RMS line profile on a velocity scale, and the error bar illustrates 
the FWHM instrumental broadening of 3 1 3 km s"' . The upper panel shows 
the value of Tccn measured in six velocity bins, and the horizontal error bars 
denote the size of each velocity bin. The solid and dashed lines mark the 
value of Teen for the entire H/3 line and its lo" uncertainty range. 



fer function (i.e., the distribution o f lag time as a func tion of 
velocity) might be obtainable (e.g.. lBentz et alj|2010bl) . 

6.3. Reverberation Lags of Ha, H^, and Hd 

We measured cross-correlation lags for the other Balmer 
lines following the same methods that were used for H/3, in- 
cluding the linear detrending. The Balmer-line reverberation 
lags follow the same pattern found previously for other AGNs, 
where the Ha line has the longest lag and higher-order Balmer 
lines have p rogressively shorter lags (e.g., Kaspi et al. 200(| 
iBentz et an i2010a). The lag of HJ is too small to be resolved 
significantly, and our measurement of Teen indicat es a la up- 
per lirn it of Teen < 2.6 days. As discussed by Bentz et alj 
( l2010al) . this trend in Balmer-line lags is the result of the dif- 
ferences in optical depths for the Balmer lines in the BLR, 
and the same behavior is seen in BLR photoionization mod- 
els jKorista & G oad 2004). For example, the Ha line has 
the highest total optical depth through the BLR, and the ob- 
served Ha emission originates primarily from larger radii in 
the BLR where the optical depth is smaller. Using the Teen 
measurements, we find a ratio of Balmer-line reverberation 
lags of r(Ha):r(H/3):r (H7) = 1.31 :1 .00:0.87, which is within 
the range observed bv IBentz et al.l (i2010a) for AGNs having 
similar emission-line lags of a few days. 

6.4. H(3 Line Width 

Virial estimates of black hole masses based on 
reverberation-mapping data employ the width of broad 
emission lines to give the typical velocity for BLR gas. Fol- 
lowing procedures often employed for reverberation-mapping 
data (e.g., Peterson et al. 2004), we measured the FWHM 
and the line dispersion (or second moment) of the H/3 line 
from both the mean and the RMS spectra. In both cases, a 
linear continuum was fitted to regions on either side of the 
H/3 line, and then subtracted off. The FWHM was measured 
directly from the continuum-subtracted Une profile, and the 
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line dispersion criine (in velocity units) was calculated as 

2 



line 



-A; 



i) h 



(3) 



where 5,- is the flux density at wavelength bin A,-, and Ao is the 
flux-weighted centroid wavelength of the line profile. 

While the narrow component of H/3 essentially vanishes in 
the RMS spectrum, it must be removed from the mean spec- 
trum before the broad-line width can be measured. The nar- 
row [O III] A4959 emission line has a nearly Gaussian profile 
with FWHM = 6.6 A, and we used this as a model to sub- 
tract off the narrow component of Hf3. To determine the final 
values and uncertainties for the lin e widths, we followed the 
Monte Carlo meth od described by Peterson et al.l (12004 1) and 
iBentz et al] (120091) . creating randomly resampled realizations 
of the mean and RMS spectra. To construct each realization, 
we randomly selected 54 spectra from our dataset, without re- 
gard to whether an individual spectrum was akeady selected, 
and then determined the mean and RMS spectra for that ran- 
domly sampled set. For a total of 1000 realizations, we mea- 
sured FWHM and criine for H/3, from both the mean and RMS 
spectra. In each realization, the borders of the continuum- 
fitting windows were shifted randomly over a range of 10 A in 
order to account for the line-width uncertainty resulting from 
the specific choice of continuum regions. The final value and 
uncertainty for each line-width measurement (FWHM and 
(Tiine) are taken to be the median value and standard devia- 
tion determined from the full ensemble of simulations. From 
the RMS spectrum, we measure FWHM = 2280 ±65 km s"' 
and uiine = 1600 ±47 km s"' for H/3. The mean spectrum 
gives FWHM = 3360 ±72 and crii„e = 1645 ± 16 km s"' for 
the broad H/3 line. After removal of the narrow component 
from the mean spectrum, the broad H/3 line centroid in the 
observed frame is Aq = 4993.2 A. The redshift inferred from 
the broad H/3 line is therefore 0.0271, and there is no signif- 
icant velocity offset between the narrow [O III] line and the 
broad H/3 emission. 

The observed line widths are affected by instrumental 
broadening in the 4" wide slit, but only by a very small 
amount. The instrumental line width adds in quadrature to 

the observed line width (as in cr^b,e,-ved = ^^Lrinsic + '^instiumental' 
or similarly for the FWHM values). We can obtain an up- 
per limit to the instrumental broadening by measuring the 
widths of arc-lamp lines observed through the 4" slit, since 
the arc lamps uniformly illuminate the entire slit width while 
the AGN does not (at least under typical seeing conditions). 
From the Cd I line at 5086 A, we find an upper limit to the in- 
strumental broadening of FWHMinst < 7.7 A, corresponding 
to an instrumental velocity width of FWHMinst < 454 km s"' 
or an instrumental dispersion of crinst < 193 km s"'. This in- 
strumental broadening is small enough that the maximum pos- 
sible contribution of instrumental broadening to the observed 
(Tiine of the broad H/3 line is only 0.9% of the observed line 
dispersion. 

To estimate the likely value of the instrume ntal disper- 
sion cTinst, we followed the method described by i Bentz et aTl 
(|2009), using an observation of [O III] through a narrower slit 
to determine the intrinsic width of the line. The initial test 
spectrum of Zw 229-015, taken on 2010 April 18, was ob- 
served through a 2" slit. In that exposure, the [O III] A5007 
fine has a raw FWHM of 5.44 A, and the 5086 A arc Une 
has FWHM = 4.18 A. Assuming that the AGN nearly fiUs 



the slit in the 2" observation (roughly the size of the seeing 
disk), this implies an intrinsic line width of FWHM = 203 
km s"' for [O III]. The instrumental broadening for the 4" 
slit can then be estimated as the difference in quadrature be- 
tween the observed width of [O III] in the 4" slit (FWHM = 
373 km s"') and the intrinsic [O III] width measured from the 
2" observation. This implies an instrumental contribution of 
FWHMinst = 313 km s-i or crinst = 133 km s"' for the 4" sHt, 
consistent with the upper limits derived previously. 

The corresponding corrections to the measured broad H/3 
line width for this instrumental broadening are very small. 
From the RMS spectrum, the corrected H/3 widths are FWHM 
= 2260 ±65 km s"' and crnne = 1590 ±47 km s"', and in 
the mean spectrum the corrected H/3 widths are FWHM = 
3350 ± 72 and crune = 1640 ± 16 km s"' . 

7. THE MASS OF THE BLACK HOLE 

Measurements of reverberation lag and broad-line width 
can be combined to give an estimate of black hole mass, under 
the assumption that the line widt h is pri marily due to virial 
motion of broad-line gas (e.g., iGaskell 1988; Kas pi et alj 
l2000t [Peterson et al.ll2004l) . From the virial theorem, we have 



MsH=f 



(cT)(AVf 



(4) 



where t is the reverberation lag and (cr) gives the mean ra- 
dius of the BLR, AV is some measure of the broad-line width 
in velocity units (typically either FWHM or the line disper- 
sion cr), and / is a scaling factor which is dependent on the 
geometry and kinematics of the BLR as well as our viewing 
angle. Since these parameters of the BLR are unknown, it 
has become customary to use a single value for / that should 
ideally represent the mean value of the virial normalization 
factor for an appropriate ensemble of AGNs. The adopted 
value of / can be based either on some set of assumptions 
about BLR kinematics (e.g., Kaspi et al. 2000), or by find- 
ing the normalization that brings reverberation-mapped AGNs 
into best agre ement with the local Mbh ~ o" relatio n of quies- 
cent galaxies dOnken et alJl2004tlWoo et al.ll2010l) . 

For consistency with the majority of recent reverberation 
work, we use criine(H/3) measured from the RMS line profile 
as the measure of AV , and cT(;en(H/3) as the measure of BLR 
size; these paramete rs have been shown t o yield the most ro- 
bust mass estimates dPeterson et al.ll2004h . For the virial nor- 
malization factor, we use /= 5.25, derived from the full avail- 
able sample of reverberatio n-mapped AGNs having measured 
stellar velocity dispersions ("W oo et a"ni2010l) . Combining the 
rest-frame value Teen = 3.86^Q gg with crune = 1590±47 km s"\ 

we obtain a "virial product" of (cr)( Ay )VG= 1.91![]i^ x lO*" 
Mq. Applying the virial normalization factor / = 5.25, this 
gives an estimated black hole mass of Mbh = 1-00-024 ^ ^0^ 

Mq. 

We note that the quoted uncertainties on Mbh above only 
include the propagated errors on the lag and line-width mea- 
surements, but not the (unknown) systematic uncertainty on 
the applied value of /. From examination of the scatter of 
reverberation-mapped AGNs about the best-fitting Mbh - o" 
relation, the typical uncertainty in the value of / applied to 
any indivi dual AGN is probably at the level of a factor of 
- 3 (e.g., Onken et al.ll2004l: IWoo et alj|2010 ). so the actual 
error in the estimated value of Mbh is dominated by the un- 
certainty in / and not by the measurement errors in the H/3 lag 
or width. A long-term goal of reverberation mapping work is 
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to use velocity-resolved measurements to directly constrain 
the BLR geometry and kinematics so that accurate black hole 
masses can be derived for individual AGNs, but such analysis 
is beyond the scope of this work. 

We can estimate the Eddington ratio L/L-eaa for Zw 229- 
015 using the V-band continuum flux and a bolometric cor- 
rection. Without high-resolution optical images or an accurate 
model of the stellar population, we cannot directly determine 
the amount of starlight contamination to the V-band magni- 
tudes or spectroscopic continuum fluxes, but a rough approx- 
imation will suffice to estimate L/i^dd- According to the ex- 
tinction map of lSchiegel et ai.l(ll99a) . the Galactic foreground 
extinction toward Zw 229-015 is Ay = 0.24 mag. We make 
the simplifying assumption that the extinction correction and 
starlight contamination correction would roughly compensate 
for each other, and take the rest-frame 5100 A flux den- 
sity from the mean spectrum as an estimate of the intrinsic 
(extinction-free) AGN flux density. We find A/a (5 100 A) 
= 4.2 X 10"'^ erg s"' cm"^ in the mean spectrum, correspond- 
ing to ALa(5100 A) = 7.1 X 10"*^ erg s"' for a luminosity dis- 
tance of 1 19 Mpc (assuming WMAP-7 cosmological param- 
eters with Hq = 1\ km s~' Mpc"'; Larson et al. 2011). Ap- 
plying the same bolometric correction used by 'Peterson et al.' 
(12004) of Lboi = 9 X ALa(5100 A), we obtain an estimated 
Lboi = 6.4 X lO'^-' erg s"', and L/LEdd ~ 0.05. This value 
is within the typ ical range found for nearby reverberation- 
mapped Seyferts (iPeterson et al.ll2004l) . 

The possible influence of radiation pressure on black hole 
masses derived from reverberation mapping has been the sub- 
ject of much recent discussion (Marconi et al. 2008; Netzer 
l2009tlMarconi et al.l2009HNetzer & Marziani20 10). The im- 
pact of radiation pressure on the motion of BLR clouds is 
likely to be strongest for AGNs radiating at high luminosity, 
but Zw 229-015 has a rather modest Eddington ratio and any 
corrections for radiation pressure should be relatively minor 
iMarconi et al.j (12008) proposed a modification to the virial 
equation, adding a term proportional to the AGN continuum 
luminosity, as in 



(crXAy)^ 
Mbh =/ ^ + i 



is 100 



1044 erg s- 



(5) 



where Lgioo = ALa at 5100 A, and where the numerical value 
of / differs from the value used to normalize the "stan- 
dard" virial equation (Equation |4]i. Using the sample of 
reverberation-mapped AGNs having measured stellar veloc- 
ity dispersions, they estimated best-fit values of / = 3.1 and 
\ogg = 7.6. Applying this formalism to Zw 229-015, we ob- 
tain an estimate of Mbh ~ 8.8 x 10^ M© forZw 229-015. This 
value is lower than the pure virial estimate for Mbh, du e to the 
combination of the smaller / value derived by Marco nTet al.l 
([2008) and the relatively low luminosity of Zw 229-015, but 
it differs from the pure virial estimate by only 12%, or less 
than Icr. From examination of the orbits of BLR clouds under 
the influence of radiation pressure, Netzer & Marzianil(t2010h 
proposed an alternative variant of the virial equation as a new 
mass estimator, dependent on the AGN continuum luminosity 
and broad-line FWHM. Applying their mass estimator (their 
Equation 18) and using FWHM(H/3) measured from the mean 
spectrum, we obtain Mbh w 1.25 x 10^ Mq for Zw 229-015. 
In either case, the corrections to Mbh for radiation-pressure 
effects are not very large, and they are far smaller than the 
overall uncertainty in the value of the virial normalization fac- 



tor / for an individual AGN. 

8. DISCUSSION AND CONCLUSIONS 

This dataset was obtained using a somewhat unconven- 
tional strategy, in which the spectroscopic observations were 
done on a nightly basis, but only during the dark runs. We 
find that this method can deliver high-quality reverberation 
results, provided that a well-sampled photometric light curve 
(not having large gaps) is also obtained. Since there are many 
facilities that can deliver nightly, queue-scheduled imaging 
observations without monthly gaps, this is now a straightfor- 
ward approach to measurement of long-term light curves and 
reverberation lags for additional AGNs. In terms of the mea- 
surement accuracy for the H/3 lag, our results for Zw 229-015 
are similar in quality to what was obtained for several objects 
with similar lag times of a few days during the 2008 Lick 
AGN Monitoring Project (Bentz et ai.i20 09). a program com- 
prising 64 mostly consecutive nights of spectroscopic obser- 
vations along with nightly photometry. 

Our results demonstrate that Zw 229-015 will be an excel- 
lent target for a variety of future investigations. Its strong 
and rapid variability is noteworthy in comparison to many 
other obje cts havin g similar black hole masses and luminosi- 
ties (e.g.. lBentz et ali2 009). and it is very fortunate that such a 
nearby and highly variable AGN happens to be in the Kepler 
field of view. While the Kepler light curve for Zw 229-015 
is not yet available, we anticipate that it will enable funda- 
mentally new investigations into the properties of AGN op- 
tical continuum variability. There are only a few other low- 
redshift Seyfert galaxies being monitored by Kepler, and al- 
though their variability properties are currently unknown, it 
seems likely that Zw 229-015 will be the most important ob- 
ject of this small sample. As such, it is crucial that Zw 229- 
015 continue to be monitored by Kepler for the longest possi- 
ble duration during the mission. 

Additional observations of Zw 229-015 can be used to mea- 
sure the galaxy's bulge properties, to place it on the cor- 
relations between black hole mass and bulge luminosity or 
mass, and to add it to the sample used to calibrate the rela- 
tionship between BLR radius and continuum luminosity fol- 
lowing the methods described by Bentz et al.l (12006). Simi- 
larly, a measurement of the bulge velocity dispersion for Zw 
229-015 can be used to add this galaxy to the Mbh-o" re- 
lation for re verberation-mapped AGNs (iNelson et al.l 120041 
lOnken et al . 2004; Woo et al. 2010). Our wide-slit observa- 
tions are not well suited to carrying out measurements of stel- 
lar kinematics, but stellar absorption features including the 
Ca II triplet absorption lines are easily visible in the Lick spec- 
tra, so it should be straightforward to measure the bulge ve- 
locity dispersion with suitable new observations. Other ob- 
serving programs for Zw 229-015 are currently in progress, 
including near-infrared monitoring with Spitzer to measure 
the dust reverberation timescale. We will also revisit our 
Hf3 reverberation measurement when the Kepler light curve 
becomes available, since the sampling and precision of the 
Kepler data may lead to an improved determination of the 
emission-line lag. 

To summarize our conclusions: we have measured a rest- 
frame reverberation lag of Teen = 3.86;!;q days for the H/3 line 
in Zw 229-015. Combining this result with the broad H/3 line 
width, we obtain a virial estimate of '--^ 10^ Mq for the black 
hole in this galaxy, under the assumption of a virial normal- 
ization factor of / = 5.25. The mass estimate does not change 
significantly if we adopt the method of .Marconi et al.. (.20081) 
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to obtain a mass estimate corrected for the dynamical influ- 
ence of radiation pressure in the BLR. The Kepler target Zw 
229-015 is destined to be one of the most important AGNs for 
investigations of the origin and characteristics of optical con- 
tinuum variability. As recent observations have begun to re- 
veal connections between optical continuum variability char- 
acteristics and black hole masses in AGNs, this estimate of 
Mbh will be a valuable asset for future interpretation of the 
Kepler light curve for Zw 229-015. 
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